Performance of DPSK Signals with Quadratic 

Phase Noise 

Keang-Po Ho* 
February 2, 2008 

Abstract 

Nonlinear phase noise induced by the interaction of fiber Kerr effect 
and amplifier noises is a quadratic function of the electric field. When 
the dependence between the additive Gaussian noise and the quadratic 
phase noise is taking into account, the error probability for differential 
phase-shift keying (DPSK) signals is derived analytically. Depending 
on the number of fiber spans, the signal-to-noise ratio (SNR) penalty 
is increased by up to 0.23 dB due to the dependence between the Gaus- 
sian noise and the quadratic phase noise. 

Keywords: phase modulation, error probability, fiber Kerr effects, non- 
linear phase noise 

1 Introduction 

Other than the projection of additive Gaussian noise to the phase, phase 
noises from other sources can be considered as multiplicative noise that adds 
directly to the phase of the received signal. When the local oscillator is not 
locked perfectly into the signal, the noisy reference gives additive phase 
noise [1, 2]. Laser phase noise degrades coherent optical communication 
systems [3-5]. Those types of extra additive phase noise that add directly 
into the signal phase are independent of the additive Gaussian noise. In 
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this paper, the additive phase noise is quadratic function of the electric 
field. When the electric field is contaminated with additive Gaussian noise, 
although the quadratic phase noise is uncorrelated with the linear phase 
noise, both non-Gaussian distributed, the phase noise weakly depends on 
the additive Gaussian noise. 

Differential phase-shift keying (DPSK) signals [6-16] have received re- 
newed attention recently for long-haul or spectrally efficiency lightwave trans- 
mission systems. When optical amplifiers are used periodically to com- 
pensate the fiber loss, the interaction of optical amplifier noise and fiber 
Kerr effect induced nonlinear phase noise, often called Gordon-Mollenauer 
effect [17], or more precisely, nonlinear phase noise induced by self-phase 
modulation. Added directly into the signal phase, Gordon-Mollenauer ef- 
fect is a quadratic function of the electric field and degrades DPSK sig- 
nal [11,14,17-23]. 

Previous studies found the variance or the corresponding Q-factor of 
the quadratic phase noise [11, 17, 24-27] or the spectral broadening of the 
signal [14, 18, 28]. Recently, quadratic phase noise is found to be non- 
Gaussian distributed both experimentally [20] and theoretically [29,30]. As 
non-Gaussian random variable, neither the variance nor Q-factor is suffi- 
cient to completely characterize the phase noise. The probability density of 
quadratic phase noise is found in [30] and used in [23] to evaluate the er- 
ror probability of DPSK signal by assuming that quadratic phase noise and 
Gaussian noise are independent of each other. However, as shown in the sim- 
ulation of [22,23], the dependence between Gaussian noise with quadratic 
phase noise increases the error probability. 

Using the distributed assumption of infinite number of fiber spans, the 
joint statistics of nonlinear phase noise and Gaussian noise is derived analyt- 
ically by [19,21,31]. The characteristic function of nonlinear phase noise be- 
comes a very simple expression with the distributed assumption [29] . The er- 
ror probability of DPSK signal has been derived with [22] and without [21,32] 
the assumption that nonlinear phase noise is independent of the Gaussian 
noise. Based on the distributed assumption, it is found that the dependence 
between linear and nonlinear phase noise increases both the error probability 
and SNR penalty [21,32]. 

The distributed assumption is very accurate when the number of fiber 
spans is larger than 32 [21, 29]. For a typical fiber span length of 80 km, 
a fiber link of 32 spans has a total distance of over 2500 km. Most ter- 
restrial fiber systems have an overall distance of less than 1000 km, the 
distributed assumption needs to be verified for small number of fiber spans. 
Recently, DPSK signals have been used in systems with small number of 
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fiber spans [16,33,34]. Of course, the independence assumption can be used 
for either small [23] or large [22] number of fiber spans. However, as shown 
in [21,32], the independence assumption of [22,23] underestimates both the 
error probability and the required SNR, contradicting to the principles of 
conservative system design. 

In this paper, taking into account the dependence between the quadratic 
phase noise and Gaussian noise, the error probability of DPSK signal is 
derived for finite number of fiber spans, to our knowledge, the first time. 
Comparing with the independence approximation of [23], the dependence 
between the quadratic phase noise and Gaussian noise increases the error 
probability of the system. 

In the remaining parts of this paper, Sec. |2] gives the model of the 
quadratic phase noise, mostly follows the approaches of [30]; Sec. El derives 
the joint statistics of the additive Gaussian noise and the quadratic phase 
noise; Using the joint statistics, Sec. 0] gives the exact error probability of 
DPSK signals with quadratic phase noise, taking into account the depen- 
dence between the additive Gaussian noise and quadratic phase noise; Sec. 
El calculates the error probability and the SNR penalty of DPSK signals, 
and compared with the independence approximation of [23]; Sec. El is the 
conclusion of the paper. 

2 Quadratic Nonlinear Phase Noise 

For an TV-span systems, for simplicity and without loss of generality, the 
overall quadratic phase noise is [17,25,26,30] 

$nl = \E + ni| 2 + \E + n% + n 2 | 2 

+ --- + \E + n l + ---+n N \ 2 , (1) 

where E$ = (A, 0) is a two-dimensional vector as the baseband represen- 
tation of the transmitted electric field, n^,A; = 1,... ,N, are independent 
identically distributed (i.i.d.) zero-mean circular Gaussian random complex 
numbers as the optical amplifier noise introduced into the system at the k th 
fiber span. Both electric field of Eq and amplifier noises of ftk in Q can also 
be represented as complex number. The variance of is -EUn^l 2 } = 2<Tq, 
k = 1, ... ,N, where cr 2 . is the noise variance per span per dimension. In 
<^Q) , the constant factor of the product of fiber nonlinear coefficient and the 
effective nonlinear length per span, 7-L e ff> is ignored for simplicity. Without 
affected the SNR, both signal and noise in (£Q) can be scaled by the same 
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ratio for different mean nonlinear phase shift of <<3?nl>= iV^4 2 + ./V(./V+l)<7o 
except the case without quadratic phase noise of <<3?nl>= 0. After the scal- 
ing, the mean nonlinear phase shift is approximately equal to the product of 
number of fiber spans and the launched power per span, especially for the 
usual case of large SNR with small noise. 

In the linear regime, ignoring the fiber loss of the last span and the 
amplifier gain required to compensate it, the signal received after N spans 
is 



E 



N 



E + ni + n2-\ h un 



with a power of P N = \E N \ 2 and SNR of p s = A 2 /(2Na%). In Eqs. (JTJ) and 
PJ). the configuration of each fiber spans is assumed to be identical with the 
same length and launched power. 

In [30], using the method of [35,36], the characteristic function of the 
quadratic phase noise Q is found to be 



N 



k=1 - 2jv<To x k 



cxp 



1 



(3) 



where w = (N, N - 1, . . . , 2, 1) T , X k , v k , k = 1, 2, . . . ,N are the eigenvalues 
and eigenvectors of the covariance matrix C, respectively. The covariance 
matrix C = M T M with 
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The characteristic function of (j3J) is used to find the error probability of 
a DPSK signal in [23] based on the assumption that the quadratic phase 
noise of Q is independent of the received electric field of ©• 



3 Joint Statistics of Gaussian Noise and Quadratic 
Phase Noise 

To find the dependence between the quadratic phase noise and the received 
electric field, the joint characteristic function of 
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**ni,,4v^'^ = E { ex PO'^NL • E N } (5) 

will be derived here with ^nl and En given by and (J2J, respectively. 
Similar to [21,30], with u = (wi,^) and i?7v = (ei,e2), we obtain 



jvNA 2 + jwi A + 2juAw T x + jioiivjx + jvzFCx, 



where 991 is given by 



(6) 



pi = |A + xi| 2 + \A + Xl + x 2 \ 2 H + |A + si H + xtvI 2 , (7) 

with nj = (xi, yi), i = 1, . . . , N, wj = (1, 1, . . . , 1) T , 

juiei = jwi(A + x 1 +x 2 -\ h x N ) = jujiA + juiwjx, 

and x = (x\,X2, ■ ■ ■ ,xn) T ■ 

Similar to [30], using the iV-dimensional Gaussian probability density 
function (p.d.f.) of (27R7q)~"2" exp ^— x T x/2o" 2 ^ for x, we obtain 



jvNA 2 +ju 1 A 



(2vr^)- 

- / exp 2jvAw T x + jiviwjx — x T Tx dx. 



(8) 



or 



^X.eiCf, ^i)=e 



(2fj5)~-det[r]-2 



x exp 



vAw + -iviwi J T 1 I zMu; + -cjiioj 



0) 



where F = T/(2ctq) — juC and T is an N x iV identity matrix. 
Similarly, using A = in ©, we get 
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tf v3>e3 (l>, W2) 



cxp 



(2cj2)-det[r]2 



where 



¥>2 = 2/i + + 2/2 1 2 H h 1 2/1 H h 2/Af| 2 - 



The joint characteristic function of 



becomes 



** NL> ^( l/ ^) = **NL( I/ )exp 



where 



(10) 

(11) 
(12) 

(13) 



*$nl(^) 



cxp 



jvNA 2 - v 2 A 2 w T T- 1 w 



(2a 2 )^det[r] 



0; 



(14) 

(15) 
(16) 



Based on the eigenvalues and eigenvectors of the covariance matrix C, 
the characteristic function of fy$ NL (i>) becomes that of ©, and 



m N (u) = A + 2jua 2 Aj2 Y 

A y ivkW)^wi)/Xk (m 
fc tt l-2jVa 2 A fc ' 1 ' 

^ = ao £l-2^A fc - (18) 

The characteristic function of ()13j) is similar to the corresponding char- 
acteristic function with the distributed assumption [21]. If the number of 
spans N approaches infinite, the characteristic function should converge to 
that of [21]. 
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Based on (j!3p. we obtain 



^ K a I = «p . (19) 

with £j, = (mjv(^),0), and denotes inverse Fourier transform with 

respect to uj. The partial characteristic function and p.d.f. of (|19|) is similar 
to a two-dimensional Gaussian p.d.f. with mean of (7Bat(z/),0) and variance 
of a 2 N (v). With the dependence on the quadratic phase noise, the variance of 
(T%(y) and the mean of mjv^) are both complex numbers depending on the 
"angular frequency" of v. The marginal p.d.f. of the received electric field 
Epf is a two-dimensional Gaussian distribution with variance of a 2 N {v)\ u= Q = 
Na 2 . and mean of w,n{i / )\v=q = A. 

With normalization, the corresponding joint characteristic of (|19|) in [21] 

has 

o 1 tan(A/jzv) 

a oo\ u ) = o 7^= — and m oo(y) = sec(Vj^)VPs ( 20 ) 

when — ► oo. Based on joint statistics of ()19|) . similar to that of [21,32,37], 
the exact error probability of DPSK signal can be derived analytically, even 
for case with linearly compensated nonlinear phase noise [23-25,27,38]. As 
shown in [21], the optimal compensation curve of [26,27] can also be derived 
using {Hi- 



4 Exact Error Probability 

With nonlinear phase noise, assuming zero transmitted phase, the overall 
received phase is 

$ r = e n - $ NL (21) 

where n is the phase of En (J2J). The received phase is confined to the 
range of [— ir, +7r). The p.d.f. of the received phase is a periodic function 
with a period of 2tt. If the characteristic function of the received phase is 
^ / $ r (^), the p.d.f. of the received phase has a Fourier series expansion of 

P*A0) = 7T + - E »{** r (m)exp(-ime)} , (22) 

m=l 

where denotes the real part of a complex number. In 1)221) . we use the 
conjugate symmetry property of ^$ r (— v) = (y). 
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In order to derive the Fourier coefficient of *$>q> r (m), we need the joint 
characteristic function of n and <I>nl at integer "angular frequency" of 
v = m. Based on 1)19(1 . using the same method as [19,21,32,39], we obtain 



X <j i ro-1 



2 

m > 0, 



+ J ro-l 
2 



7(l/) 



(23) 



where 7(1/) = imjvMV^ArW ^ s the complex-valued frequency dependence 
SNR parameter. When v = 0, it is obvious that 7(^)^=0 = p s - 

From (|2*T|) . the Fourier coefficient in ((2*2l is ^$ r (m) = ^$ NLi e n (m, m). 
For DPSK signal, the differential received phase is A$ r = <& r (i) — $ r (t — T) 
in which the p.d.f.'s of 3> r (i) and $ r (t — T) are the same as that of (1221) . 
The p.d.f. of the differential received phase is the same as 1)22(1 with Fourier 
coefficient equal to |^$ r (m)| 2 , i.e., 



^ ^ +00 

P^M = 1T + -H \^*r(™)\ 2 cos(m0). 



(24) 



m=l 



Similar to the procedure of [2,3,21,23,32,39-41], the error probability 
becomes 



1 



+00 



9 TT ^ 



-1) 



|¥* r (2fc + l) 



(25) 



or 



Pe 



k=0 



-l) k \r k e rk \ 
2k + 1 



r k 



4hf + 4+1 tt 



x|^ NL (2fc + l)| 2 , 



(26) 



where 



m 2 v (2fc + 1) 



(27) 



2(7^(2^ + 1) 

analogous to the "angular frequency" depending SNR as the ratio of complex 
power of ^m 2 N {v) to the noise variance of a 2 N (v). 
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Figure 1: The error probability of DPSK signal as a function of SNR for 
N =1, 2, 4, 8, 32, and infinite number of fiber spans and mean nonlinear 
phase shift of <<3?nl>= 0.5 rad. 



The error probability expression of l)26|) is almost the same as that in 
[21, 32] but with a different parameter of (|27ft . The error probability of 
(|26j) is also similar to the cases when additive phase noise is independent to 
Gaussian noise [2,3,23,40,41]. The frequency depending SNR is originated 
from the dependence between the additional phase noise and the Gaussian 
noise [19,21,32,37]. 



5 Numerical Results 

For DPSK signals with quadratic phase noise, Figure ^ shows the exact 
error probability as a function of SNR p s for mean nonlinear phase shift of 
<^nl>= 0.5 rad. Figure El shows the SNR penalty for an error probability 
of 10 -9 as a function of mean nonlinear phase shift <<3?nl>- The SNR 
penalty is defined as the additional required SNR to achieve the same error 
probability of 10~ 9 . Both Figs. Q and [21 are calculated using (f2*B|) and the 
independence approximation of [23]. The independence approximation of 
[23] underestimates both the error probability and SNR penalty of a DPSK 
signal with quadratic phase noise of Both Figs. HandElalso include the 
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Mean Nonlinear Phase Shift (rad) 
Figure 2: The SNR penalty vs. mean nonlinear phase shift <<I>nl>- 

exact and approximated error probability for N = oo that are the distributed 
model from [32] and [22], respectively. The distributed model is applicable 
when the number of fiber spans is larger than 32. In Fig.^ without quadratic 
phase noise of <<1?nl>= 0, the error probability is p e = exp(— p s )/2 [42]. 
The required SNR for systems without nonlinear phase noise of <<3?nl>= 
is p s = 20 (13 dB) for an error probability of 10 -9 . 

From Figs. ^and[2J for the same mean nonlinear phase shift of <<3?nl>) 
the SNR penalty is larger for smaller number of fiber spans. When the mean 
nonlinear phase shift is <<3?nl>= 0.56 rad, the SNR penalty is about 1 dB 
with large number {N > 32) of fiber spans but up to 3-dB SNR penalty 
for small number {N = 1,2) of fiber spans. For 1-dB SNR penalty, the 
mean nonlinear phase shift is also reduced from 0.56 to 0.35 rad with small 
number of fiber spans. 

In [17], the optimal operating point is defined when the variance of 
quadratic phase noise is approximately equal to the variance of the phase 
of Gaussian noise. In [22,23], the optimal operating is calculated rigorously 
at the operation condition in which the increase of launched power does not 
improve the system performance. The optimal operating point is reduced 
from 0.97 to 0.55 rad with the decrease of the number of fiber spans. 

When the exact error probability is compared with the independence 
approximation of [23]. The independence approximation is closer to the 
exact error probability for small number of fiber spans. In all cases, the 



10 



independence assumption of [22, 23] underestimates the error probability 
of the system, contradicting to the conservative principle of system design. 
The dependence between linear and nonlinear phase noise increases the SNR 
penalty up to 0.23 dB. 

Prom the SNR penalty of Fig. |5J if a prior penalty of about 0.23 dB is 
added into the system, the independence assumption of [23] can be used to 
provide a conservative system design guideline. 

6 Conclusion 

For a system with small number of fiber spans, the exact error probability 
of a DPSK signal with quadratic phase noise is derived analytically the first 
time when the dependence between linear and nonlinear phase noise is taking 
into account. For the same mean nonlinear phase shift, the error probability 
increases for small number of fiber spans. The dependence between linear 
and nonlinear phase noises increases the error probability for DPSK signals. 
Depending on the number of fiber spans, the SNR penalty increases by up 
to 0.23 dB due to the dependence between Gaussian noise and the quadratic 
phase noise. 

For the same mean nonlinear phase shifts and SNR, the error probability 
of the system increases with the decrease of the number of fiber spans. As an 
example, the optimal operating point for system with large number (N > 32) 
is a mean nonlinear phase shift of about 1 rad that is reduced to about 0.55 
rad for system with small number of fiber spans (N = 1,2). 
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